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ABSTRACT 


Aperture averaging effects on the measurement of 
Scintillation and on the temporal-frequency power spectrum 
have been examined in the marine boundary surface layer 
using a He-Ne (6328 A) laser beam propagating over various 
path lengths and under varying turbulence conditions. Two 
Significant results have been obtained: (1) For scintilla- 
tion measurements, a vanishingly small receiver can best 
be approximated by using an aperture of diameter 


Dp... = (z/k)?/?. 


meee (2) The aperture averaging correction 


factor for an aperture diameter D > Dae decreases as the 
level of turbulence increases, 

Data have also been obtained supporting the theoretical 
prediction that an increase in aperture diameter reduces 


the high frequency power content of the temporal-frequency 


power spectrum. 
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I. INTRODUCTION 


Point iibartOy Lsethe term used to describe the 
atmospherically induced intensity fluctuations of a constant 
power output laser source as seen from a distant detection 
device. The precise manner in which a detector's aperture 
Size affects the measured value of the variance of the 
Mmronsity tluctuations as a function of distance from the 
source and lasing wavelength has not been satisfactorily 
established. Several aperture averaging experiments have 
been conducted under various turbulence conditions to examine 
this effect along with an analysis of the frequency power 


spectrum effects of aperture Size. 


hee ROJECT HISTORY 

Since 1971 an interdisciplinary group of professors 
principally from the Departments of Physics and Chemistry, 
Meteorology, and Mechanical Engineering have been investi- 
gating under a Navy research contract various aspects of 
laser beam propagation over the marine boundary layer. Data 
has been gathered thus far on (1) Modulation Transfer Function 
(MTF) at wavelengths of .4880 um, .6328 um, 1.06 um and 10.6 un, 
emccinttiiation at .488@ wm, .6328 wm and 10.6 pm, (3) extinc- 
tion at .6328 um, (4) beam wander and (5) beam spread. 

The stated general objectives of the project are: (1) to 
develop an understanding of the factors affecting the propaga- 


tion of optical radiation through the atmosphere, particularly 
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the marine boundary surface layer, the lowest 30 meters 
above the ocean surface; (2) to develop the capability of 
predicting optical system performance on the basis of point 
meteorological measurements. 

Previous thought had been that the sea surface might 
pmesent a smooth, isothermal layer, giving rise to little 
turbulence effect on a propagating light beam. However, 
early experiments in the project proved that this was not 


the case. 


ime Oo GANT ILLATION 
im Description 

The twinkling one sees when observing distant light 
sources is a physical realization of scintillation. This 
appears as fluctuation in the output from a small photo- 
detector illuminated by a distant source. 

In the general case of finite beam propagation, 
measured intensity fluctuation can be caused by any or all 
of the following: (1) source output power fluctuations, 
(2) beam wander, (3) detector motion through a plane 
perpendicular to the optical axis, (4) a temporal wavefront 
breakup within the beam caused by small-scale spatial and 
temporal variations of the dielectric constant (or index of 
refraction) of the lower atmosphere. Of these, the most 
critical for experimentation or information transmission are 
the turbulence effects of scintillation, beam wander and 
beam spread, since they cannot be overcome by system 


refinement. 
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2. The Effects of Turbulence on a Propagating Wavefront: 
ysical an eoretica odels 


a. Kolmogorov Model of Turbulence 

Currently the only well-developed model of 
atmospheric turbulence is that of Kolmogorov. It is basically 
a velocity turbulence model that assumes homogeneous, iso- 
tropic turbulence. The turbulence is assumed to be composed 
of small turbulent eddies, or turbules, each of characteris- 
tic temperature and density (hence index of refraction), as 
illustrated in Figure 1. In this homogeneous turbulence 
model, often referred to as "frozen-in turbulence," fluctua- 
tions in the dielectric constant over a vanishingly small 
volume element are due to atmospheric temperature variations 
produced by the turbules being swept past the reference 
point by the mean wind. 

Fete m(t,t) be the reiractive index at point r 


in Space at time t. Then 


G 
=f atr,e)de + my (r,t', 2) 
0 


(1) Imig) 


<n(r,t)>_ + n,(r,t"',7) 


where n,(r,t',7) is the deviation from the average value at 
point r and time t. 
The refractive index structure function gives 


the mean square variation in the refractive index between 


de? 





BvorpOlntS in space: 
2 
(2) DC eem in) = n,(25)) > 
where oS Ir, - r,| 


The cross-correlation function of refractive 


index fluctuations is defined as 


cL 
(3) By(e) = Ef ny (eyet ay (eye! )ae! 
0 


<n, (4, Jn, (r2)> 


Under the assumption that homogeneous turbulence is present 
and the appropriate statistical averages for BL Ce) and 
D_ (0) are Stationary in time and locally constant in space, 


one finds that 


(4a) D.(e) = 2{B (0) - B(e)] 
(4b) B_(e) = 5{D,(#) - D,(0)] 


As 9 becomes smaller, correlation increases while the 
structure function decreases. 

According to the Kolmogorov theory of turbulence, 
for wae =< Li, that is, within the "Inertial Subrange" of 


Readies, 
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(S) D,(e) = Cy” p ; 


and the corresponding temperature structure function for 


meteorological measurements is 
(6) Dr (0) = Cp? p2/5 
- ale 


de and Lo are the inner and outer scales of turbulence 


N and Cor are the 


refractive index structure constant and temperature struc- 


in the Inertial Subrange (Figure 2), and C 


Mimenconstant, respectively. If % (or L) represents the 
diameter of a turbulent eddy, de is the scale size below 
which viscous effects cause the conversion of turbulent 

energy into heat, and subsequent destruction of the eddy 


ensues. L the outer scale size, is the minimum scale 


6? 
at which energy enters the turbulent system from wind shear 
and heat flux. 

For Kolmogorov turbulence, an analytical form 
(the von Karman form) of the power spectrum (the three- 
dimensional Fourier transform of the cross-correlation 


function), which illustrates its approximate dependence on 


do and Lo 1s 


-11/6 


2) 


(7) 9% (K) = 0.033 C7 (Ke 2 5p exp [- (K2,) 7] 
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hon K >> Lo» this reduces to the more common form used by 


Tatarski, 


E . Ze ef 3 2 2 
(8) on) O7033 C. K exp [-K 1 ee Ene: 
where K = 2T 
m de 


b. Physical Models 
Hiimelatarski's Geometrical Optics Approach. The 

physical mechanism through which turbulence induces beam 
breakup is by the focusing of sections of a wavefront 
produced by the inhomogeneities in the index of refraction 
over the propagation path, Tatarski uses a plane wave 
incident upon a series of closely packed refractive spheres 
of various sizes. Using the paraxial approximation, a 


Sphere acts as a lens with focal length f = x > where n, 


1 
is the index of refraction perturbation within the sphere. 


For a = Im, Xo of the Inertial Subrange, and n, = Tone ; 

££ = 500m. For a = 1m, the lowest L, lime, tf = Ssv0 Km. 
Thus, for most applications in a turbulent medium, a detector 
will appear to be closer to the "lens" than to the focal 
plane. The fractional intensity change induced per sphere 
will be small. Over a long path, perturbation effects will 
be multiplicative and the resulting wavefront at the 


receiver will be neither coherent nor of conStant intensity. 


ineaddwerons dilitraction effects, important at distances 


He 





greater than the Rayleigh range, Z, = ka’, limit the ability 


R 
of the sphere to focus the incident wave. The amplitude 
fluctuation induced by a sphere at a distance Z from a 


detector is estimated [Ref. 1] as 


2|n,{zZ 


O 


From equation (9) it is seen that the largest fluctuations 
are produced by the smallest turbules. By the geometrical 
optics approach [Ref. 1], within the Rayleigh range, 
amplitude fluctuations increase as the cube of the distance, 
are dominated by the smallest turbules, and are independent 
of wavelength: 


ba 2 _ 
(10) fe ee” a? 7? 
O 


However, this model does not take into account diffraction 
effects, important for smaller turbules that contribute most 
toward amplitude fluctuations. 

Camelia tiracetom Effects.) Diffraction effects 
are taken into account by assuming that spheres of radius 
1/2 


gee~ (AZ") 


A , Where Z' is the distance from the center of 


the turbule to the receiver, contribute very weak, diffracted 
waves with no focusing effects at the receiver. Substituting 
(azyh/2 OTs de in equation (10), the mean square amplitude 


fluctuations become 
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. ke? SOARES : 


ge 


C 2 , 7/6 711/6 


SA, 2 
(11) <q N 


Now defining a a where Xy is) the first order jlog- 
O 


ampiitude for |X,| << 1, Tatarski has calculated the first 
order log-amplitude variance with diffraction effects con- 


Sidered for plane waves as 


(12) <x, 7 = 0.31 ,//6 eS 7i1/6 
and for spherical waves as 
(13) <x,*> = 0.13 Ki 16 a on 


both restricted to Z >> ke? Lutomirski and Yura [Ref. 1] 


State that one would expect the behavior of a beam wave to 
lie between that of the plane and spherical waves. There- 
fore, one may use equations (12) and (13) as bounds for 
weak scintillation. For inhomogeneous turbulence, equations 


(12) and (13) can be expressed in the form 


AAOD 


(14) <X4 


7 
Sesame 2 f (z-2') 2/6 Cy (Z") dz! eed 
Oo 


7 
(15) <x 17> ©) rOnGcme © f (z-z')>2/ ae Cy (21) dz" 
O 
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where Z is the distance from source to receiver and Z'! 
from turbule to receiver. For plane wave propagation, the 
fluctuation effects are weighted by the factor (ion 
which has its maximum value near the source and approaches 
zero at the receiver. The spherical wave weighting factor, 
Me-2')(--)1°/°, has itS maximum value at the midpoint of 
the path and approaches zero at both ends. 

c. Theoretical Models, The Scalar Wave Equation 

Numerous attempts have been made to solve the 

Scalar wave equation to relate the electric field E(r) 
of an optical wave passing through a turbulent medium to 
the refractive index fluctuations. Three approaches current- 
ly in use for solving this problem are: (1) Geometrical 
Optics, (2) The Born Approximation, and (3) The Rytov 
Approximation. Below is a brief resumé of the approaches, 
their assumptions, highlights, methods of approach and 
restrictions. All three methods start with the scalar wave 


equation, 


Z 2 


(16) V°E + k“n“E = 0 

The derivation of this has assumed a medium with no free 
charges or current density, with the permeability of free 
Space, that correlations between the components of the 


field polarized in different directions can be neglected, 


and that A << Lee 
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(i }@Geometricale@ptics Approach 

(a) Amplitude fluctuations are attributed 
temtne £Ocusing and defocusing of the rays by the curvature 
of the turbulent eddies along the path. 

(b) Phase fluctuations are calculated from 
the changes in velocity as a ray passes through regions 
of different refractive index. 

(c) Diffraction effects are ignored. 
t7 2 


<< 2 


(d) The major restriction is (AZ) a 


This reduces the applicability of this model to extremely 


1/2 


Short ranges. Note: When (AZ) > & diffraction effects 


0? 
become important. 
(2) The Born Approximation 
(a) This 1s a perturbation technique to 
incorporate diffraction effects into the scalar wave equation. 
teen + ny is substituted into the scalar wave equation, 


woaten 1S solved for the first order solution for E. The 


result is 


(17) Ey - [ 2k*n, (x'JES(e')GCl[z-z'|Jav' , 
V 


PibemenG(jr-r'|) is the free space Green's function. 

(b) The perturbation technique used is valid 
as long as the amplitude fluctuations remain small. 

(c) It is assumed that the field at the 
receiver is the sum of the original incident field and the 


fiéld Scattered one time from a turbule. 
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(d) It requires that (az) 1/2 > he 
(e) Restrictions: In the optical wavelength 
range, amplitude fluctuations often become very large 
compared to the average value of the amplitude. Multiple 
Scattering effects cannot be ignored (but are). References 


2, 3 and 4 contain details on multiple scattering effects. 


(3) The Rytov Approximation (Method of Smooth 
Perturbations ) 


(a) This is a perturbation technique applied 
to a transformation of the scalar wave equation using the 


exponential notation: 
e(r) = e¥') | yer) = x(r) + id(r) 


w(r), the "complex amplitude," contains both the amplitude 
and phase terms. Using this, the scalar wave equation takes 


the form of the non-linear Ricatti Equation, 
2 Z = 
(18) Vp + VweVbd + k°(1 + On) = 0 


Now applying the perturbation technique ~W = y+ Vi> where 


O 
Vy is the perturbation to the average value of Vo» one 


obtains the solution, 


a9) vr) = ty f cern) vy, +1, + Kon]E,CF)av" 
O 
V 
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ine “Rytov Approximation" is Ve °V¥, = Ivy, 17 = 0, and 


results in 


. 2 it 
(20) Yip (4) = Soe! G(r-r')én(r')E,(x')dv' , 
O 
V 


tiem tirst order complex phase.'t A higher order solution 
On Y, can be obtained by iteration. 

(b) This approach includes multiple scattering 
effects. Caveat--these effects are included in an inflexible 
way which does not depend on the turbulence or other corre- 
lated parameters [Refs. 5 and 6]. 

(c) The method requires that the perturbations 
of the log-amplitude, the real part of the complex phase 


[X, = log. (s~)], and of the phase, by» be small, i.e., 
O 7 


A|X, | << 1 and A[9,] << 1. 
(d) Restrictions: A << 2, IVY, | << ik76n], 
3 
AL 
— abe 
O 


(e) This method includes both the Geometrical 
Optics approach and Born Approximation as special cases. 

(f) Theoretical evidence indicates the 
equation derived for amplitude fluctuations is valid as 
long as <y “> elmo odturacron Condition) (Refs. 7-10]. 


3. The Form of the Probability Distribution for 
Amplitude Fluctuations 


The application of the Central Limit Theorem in 


the derivation of the log-amplitude equation leads to a 


Zt 





prediction of a log-normal distribution of the amplitude 
fluctuations, or equivalently, a normal distribution for 
log-amplitude. Some authors [Refs. 11 and 12], taking a 
different approach, argue for a Rayleigh distribution. The 
Rayleigh distribution arises when the field at the receiver 
mchne sum of the electric field components from randomly 
Scattered wavefront elements. This occurs in line-of-sight 
propagation when the effects of the turbulent medium are 
far from the receiver. In line-of-sight propagation where 
the receiver is located within the region in which turbulent 
effects are being introduced, the field at the receiver is 
the result of multiplicative effects on a single wavefront. 
Strohbehn [Ref. 2] shows that v(t), representing the ampli- 
tude and phase terms of a single wavefront, is composed of 
the sum of a large number of independent components, and 
from the Central Limit Theorem, has a normal or Gaussian 
distribution. 

For long enough paths, however, both the multiplica- 
tive effects on a Single wavefront and additive effects 
from different wavefronts may be important, and Strohbehn 
indicates that it is difficult to tell which effect will be 
predominant. Empirical results thus far support a log- 
Moumal distribution [Refs. 14, 15 and 16 and Figure 9]. 
However, extremely long-range experiments have not as yet 


been conducted. 
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C. THE SATURATION OF SCINTILLATION 
1. Description 
For plane and spherical waves (only) in homogeneous 


turbulence and for weak fluctuations, 


ee 
Zane ~ 2 ae Z 
21 ) an es 4 SX a = ora ; 
where a is the normalized variance and aoe the log- 
intensity variance. In terms of the log-intensity variance, 


equations (12) and (13) become 


7/6 ~11/6 2 
aI eee Kk Z Cy 


(22a) O (Plane Wave) 


(22b) o 0.52 k 


716 zil/6 ne (Spherical Wave). 
Empirical results [Refs. 14, 32-34] have shown that when the 


optically measured To oa 1s plotted against Tatarski's 


I? 
theoretical values of equation (22), g,7 8%), calculated 
from the temperature structure constant, the measured 
Scintillations reach their maximum value for Sek? of the 
Order of 1 to 3, then decrease slowly and tend to level off 
With increasing fe al This approximately constant value 
ee 1 See 2 1 is called the saturation effect. Yura 
[Ref. 17] states that the log-amplitude covariance function 


differs both quantitatively and qualitatively in the 


Saturation region from that under unsaturated conditions. 


Zo 





The actual value of ta 


Masemot yet been strictly established. Strohbehn [Ref. 2] 


at which saturation begins 


indicates that saturation begins for o os 0.8, while 





QI 
de Wolf [Ref. 18] and Lutomirski and Yura [Ref. 1] maintain 
Z 
mroccurs for Soy ? ik 
2. Theory 
De Wolf has defined the saturation region in terms 
ars ie as 
2 kL,” 
(23) = eee eS 
or equivalently, 
2 
2 OSES 


(24) ee kZ a. ie 
De Wolf [Ref. 18] improves upon the Rytov approximation by 
Pirnodiciwes a directional correction factor, cos[K-p(Z)], 
to allow rays being refracted by a turbule to be deflected 
by more than one Fresnel zone size off the optical axis. 
This directional factor eliminates the focusing effect of 
turbules close to the detection point and introduces serious 
interference effects and therefore amplitude changes. 

De Wolf has calculated the log-amplitude variance in 
the saturation region based on plane wave theory to be 

2 3 Dee 6 27 


: 7/3 Z 
(25) <x, °> = 0.41(K “4° 2° C,”) = ae 
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and derives the result that the irradiance distribution 
remains log-normal in the saturation regime. However, to 
verify this equation, Kn must be determined empirically, 
Which is Foe eene COmdor 

Other authors have attempted to describe the 
Saturation effect by carrying higher order terms through the 
Wave equation. The efforts have been extremely complex 


mathematically and have resulted in little success. 


D. TEMPORAL-FREQUENCY SPECTRUM ANALYSIS 
ieee Description 

Like any other time-varying signal, a fluctuating 
voltage response from the output of a photodetector possesses 
an autocorrelation function. The power spectrum of the 
Signal can be found (1) manually by computing the Fourier 
transform of the signal's autocorrelation function, or 
(2) electronically by feeding the fluctuating voltage signal 
into a Spectrum Analyzer, which performs the transform. The 
power spectrum shows the frequencies composing the signal 
and their relative magnitudes. Current theories (described 
below) indicate that the maximum frequency that will be 
present will be directly proportional to the mean wind speed 
normal to the propagation path and inversely proportional 
to the square root of the wavelength used and distance of 
Peopacation. Theories differ significantly on the effects 


of aperture diameter on the power spectrum. 


Ke. 





2. Theories 
a. Lutomirski and Yura 
It 1s assumed that all refractive index 
inhomogeneities of the turbulent medium move convectively 
With mean wind speed v everywhere along the propagation 
path. Intensity fluctuations will be dominated by the 
velocity component normal to the propagation path, Vn If 
the transverse intensity correlation length is Po» then the 
correlation period of the intensity fluctuations at a point 
is of the order of we P/V,» giving a characteristic 


frequency of 





(26) esa 


Using Lutomirski and Yura's [Ref. 1] model of Py ~ (z/k) 1/2, 
the maximum frequency cutoff for scintillation effects at a 
point is given by 


V 
OF Bo = 2 
-? © ¢am12) 272 


The normalized photodetector signal at time t 


1s given by 


I(P,,t) aes 


(28) S*(t) = oO) a = ee ae 


and from the frozen-in turbulence assumption, 
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(ea Vet. c) = <i> 
(29) S4(t +t) = f —4A Ba? 


Mieeeuctocorrelation function of S'(t) is 


= ! ! = 1 7 = 2 2 
fo) R_,(t) = <s'(t)s'(t+1)> es Gra |e ca eid ge 


The power spectrum is the Fourier transform of Ro, (t) and 


co 7 aN ence 
2 L’wy 2 
pea es 8) ae 


(31) 6,4) = “2 | c,Avfiex? ot 


<I>" w 0 Vn \ /1+x? 
D 


Aperture averaging effects on the spectrum are seen through 


Hoegiven by 


ale 


JN 
the Wy term, wy = + Where D is the aperture diameter. 


Jy is a first-order Bessel function and C, (K) 1s the two- 
dimensional power spectrum of the intensity fluctuations. 
The square bracketed term falls off rapidly for w > Wh» 
displaying the filtering action of the aperture diameter. 

When plotting a normalized power spectrum, 
o,(£), versus log f, a change in A, Z, Wis or D effects a 
shift of 6, (2) along the log f axis. 

beeeclatrord [Rer. 19] 

A theoretical solution for the temporal- 

frequency spectrum of scintillation at a vanishingly small 


aperture is given as shown for plane and spherical waves 


in Figure 3. The corresponding asymptotic forms are given 


a7, 





for plane waves as 


? 


(32) (f) = 0.44 (6..°) G4 (1 + 0,27 af/3) a9 xc 
: O 


8/3 


(33) O(£) = 1.44 (c,,-)(2) 2? , a> 1, 
. O 


and for spherical waves as 
Zee 473 
34.) CeO. 245 (Ton (gy Tt ve (OFS a t9. / ls < le 
fe oe 8) 5 
(35) DE) = 258 1 (Fon (er) Q pease ae a, 


f — -1/2 
where 2 = F- and fA = V (27 AZ) / : 


O 
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Puke OF THE PROBLEM 


A. APERTURE AVERAGING EFFECTS ON SCINTILLATION 

1. Description 

A photodetector placed at the image point of a 

lens illuminated with a plane wave responds to the square 
of the integrated amplitude over the convergent wavefronts. 
Having passed through a turbulent medium, the wavefront of 
incoming light for an ideal spherical or plane wave is no 
longer coherent nor of equal amplitude across the surface 
at the detector aperture. The random-sized turbules 
of different refractive indices along the propagation path 
have perturbed the wavefront creating a buildup of amplitude 
in some areas and a cancellation in others. A transverse 
eomrelation distance, Po; is defined as the distance in which 
intensity fluctuations remain correlated on the traveling 
wavefront at the point of observation, and in general is 
a function of path length Z and wavelength of the source. 
If a detector of diameter D is used to measure the intensity 
fluctuations at a point distant from the source, then the 
variance of the log-amplitude measured will depend on the 
value of D as compared to Pos On ol Pos the wavefront 
incident upon the detector will not be strongly correlated, 
cancellations from the amplitudes of different parts of the 
wavefront being out of phase will occur, and intensity 


fluctuations measured will be less than for a correlated 


Ze 





Moe rront, 1.e., D < Po. This observed decrease in the 
variance of the log-intensity fluctuations due to increasing 
aperture size is called the aperture averaging effect. 
2. Theories 
Tel MtOlirck! and Yura [Ret. 1] 
For small scintillations and homogeneous 
turbulence, the transverse correlation distance of the 


intensity fluctuations is given as 


(36) Ay = (2/k) 1/2 = .399 (az)1/? 
For strong scintillation, no reliable estimate of Po exists. 
If the aperture diameter D < Pos then the intensities 
measured at all points ry Within the aperture, I(r,,t), are 
correlated at time t, and the photodetector current s(t) 


1s proportional to 
x 1D¢ 
(37) s(t) = = 1(0,t) , 
1l.e., an effective point detector. As D becomes greater 
than Pos s(t) fluctuates less and for D >> Po, the average 


incident intensity becomes constant and we have 


(38) <S> = mD <I> 





The signal variance is given by 
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2 Zz Z 
(39) Be = fertty-mb4 p,d Ps 


where C, (Pp) is the covariance of irradiance. Lutomirski 


and Yura state that aperture averaging ""...is the ratio of 


a.* in (equation 39) to what the variance would have been 


if the intensity over the whole aperture varied as it did 
mea point.’' At a point, 
ee 
om 711) 
(40) 0 = Fae C, (0) 


> 
O 


The aperture averaging factor, A, is defined as: 


Z 
2 GG py ps1) 
_ Ss - 4 ee 7 Z 2 
(41) = a Sa f Oa d pd Eo) 
Ss 
O 


Reference 1 shows that if: 


i 7 D, then A = 1, i.e., there is no aperture averaging 
mae tOr ; 
~ Ph, é 
on), elena = 4 (5) : 


Po ~ D, then the value of A requires a detailed knowledge 
of Cr (Pp). C; (CP) varies with source diameter, propagation 
Slant path and meteorological conditions. For homogeneous 
turbulence and plane wave propagation, Cr (Pp) is not known 


for strong scintillation. It is recommended that 


SL 


7 — =. =. == 





(42) A = —__—__, 
Lo Ga) 
R 


be used for homogeneous turbulence when the scintillation is 
small. 
b. Fried [Ref. 21] 
Using the statistics of propagation of an 
infinite plane wave over a horizontal path, Fried derives 
equation (37), then develops an expression for the signal 


variance, 


D 


(43) 0.7 = 2n J 0K (p,D) Cy(o)do , 
O 


where K, (9, D) is the area of overlap of two circles of 
diameter D whose centers are displaced a distance p, and 
is given by 


2 
eee jie oc Leen 
(44) Ki (o,D) = ¢ 4 ee = 


Fried then derives the same aperture averaging factor as 
given in equation (41) but labels the factor 9. He states 
that the normalized log-amplitude covariance, C,(e)/C, (0), 


sel e2 


Pomamiimecion of o0(42/k) only, for the horizontal 


propagation case, and gives for an aperture diameter 


OZ 





normalization factor, 
(45) (4z/k)1/4 eS (az) 2/2 


Fried's theoretically calculated dependence of 
Q upon the normalized aperture diameter, D/ (42/k) 1/2, 1s 
shown in Figure 4 for various values of C,(0). 
¢. Tatarski [Ref. 20] 
This theory on aperture averaging effects 
gives the amplitude correlation distance as being "on the 


2 


order of" the first Fresnel zone size, (AZ) POd oe els 


model then, the aperture normalization factor is simply 





cazyi/2. Tatarski was the first to define the aperture 
2 
. = Q(R) _ Ss 
averaging factor G(R) = OCOy ° where Q(R) = ee and 
o. é 
Q(0) = =? , Which Lutomirski and Yura have labeled "A" 


and Fried "9," 
3. Published Experimental Data 
a. Homstad and Heneghan [Ref. 22] 
Measurements were made on a mesa near Boulder, 
Colorado, over a horizontal (statistically homogeneous) path 
at distances of 200-1600m using a He-Ne (6328 Ay laser, 


G(R) vs. R/(Az)1/? 


was compared with Fried's theoretical 
curve of Figure 4. G(R) is the aperture averaging factor, 
equivalent to A in equation (41). The result was that good 


agreement with theoretical calculations was observed 


a 





" ..in the region of weak fluctuations (i.e., at shorter 
path lengths)...data in the region of strong fluctuations 


...show a trend toward larger G(R) for R/ (az) 1/2 


a ee 
agreement with eee [Rets. 235 and 24]. This is consistent 
with recent theoretical predictions [Refs. 25 and 26] and 
with Kerr's experimental work which found the correlation 
radius to be longer in the region of strong fluctuations 
Meets, 25, 24, 35 and 36]. 

Homstad and Heneghan point out that a reason 
that experimental data collected by Fried [Ref. 16] did not 
support his theoretical predictions may have been that the 
apertures he used for the 8 km path length were not large 
enough to extend beyond the correlation radius and hence 
gave a relatively constant G(R) as reported by Fried. 

b. Strohbehn and Homstad [Ref. 27] 

In an experiment conducted over land, good 
agreement with Fried's model in the region of weak fluctua- 
tions is confirmed. 

c. Kerr [Ref. 23] 

States that the predicted value of the aperture 

averaging factor should vary as D ~ for diameters much 


larger than (az) 1/2 


, Where different authors have assigned 
values of 2, 7/3 or 3 to x, and that he did not observe this. 

d. Previous NPS Work [Ref. 28] 
Schroeder observed that "Significant aperture 


averaging was observed in an enclosed laboratory for 


apertures less than the first Fresnel zone size."' From 


34 
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two field experiments (both conducted at the same path 
length and with the same wavelength) Schroeder stated that 
the aperture averaging factor may be dependent on Cy and 
mee *'...a 4.8 mm aperture can approximate a vanishingly 
small receiver and should be used whenever possible.'' The 
two field experiments were conducted transmitting across 
Monterey Bay with receiver and source located on shore. 
Path length was 4.05 km, 98 percent of which was over the 
ocean surface at 4.9 meters. Aperture sizes used varied 
Peo. > CM tO 2.35 cm. By graphing "<9, (D) - Oo, , (23mm) > 
vs. Aperture Diameter D,' Schroeder's data indicated Thee 


the largest average difference in o A values was seen for 


x 


a 4.8 mm aperture. 
Peer eR TURE AVERAGING EFFECTS ON THE TEMPORAL- FREQUENCY 
POWER SPECTRUM 
1. Tatarski's Theoretical Model 
With increase in aperture size, Tatarski predicts 
a loss in the high frequency components of the normalized 


Spectrum. By the addition of the factor 





Z 
io : ane ee/y 2 
RVK 2 + An £2/y 
nh 


in the kernal of the integral for the plane wave model of 


W(f), the spectral density of intensity fluctuations, a 
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partial suppression of the normalized intensity fluctuations 


as seen by an aperture of radius R is accounted for. For 


>? sr EES correction factor is a very small number, 
which corresponds to a suppression of high frequencies. 
As the aperture size is increased, the frequency at which 
the suppression occurs decreases. The low frequency portion 
of the spectrum will remain unaffected [Ref. 20]. 

Using the frozen-in turbulence model, the relation- 
Ship between large spatial frequencies, K, and large temporal 


frequencies is given by w = Kev, . This implies that a 


receiver of diameter D will pass only temporal frequencies 


<| 


of w <r : 
Z. Published Data 
a. Hohn [Ref. 29] 

A He-Ne (63284) source was used over path lengths 
of 4.5 km and 14.5 km. Various aperture sizes from .5 cm 
to 8 cm were used. Wind speed values were obtained from 
a meteorological center near the test site. Normalized 
Spectrum Curves were deduced from the data taken and the 
results were that "No systematic correlations of U(f)... 
or fo/2 were found with changing aperture diameter D of the 
photometer."' A second result was that the theoretical 


zi/2 did not appear to hold empirically. 


rwelation f. x 
b. Schroeder's NPS Work [Ref. 28] 
An aperture averaging frequency analysis sequence 


was conducted during each of the two field experiments 
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previously mentioned. Schroeder noted that changes in wind 
velocity could produce a different spectrum if two data 
samples at different aperture sizes were not conducted 
simultaneously. General wind speed conditions were measured 
at the transmitting location and assumed to be the same 

over the propagation path. The spectrum was obtained by 
passing a recorded scintillation signal through a wave 
analyzer with a bandwidth of 10 Hz at 10 Hz intervals from 

20 to 100 Hz, then at 100 Hz intervals from 100 to 1000 Hz. 
The power was determined by integrating the output curves 
with a planimeter. Results were that for f < 90 Hz, the 
aperture averaging effect on the relative power spectrum 

was "unpredictable," while for £ > 90 Hz, the smaller 
apertures show a larger contribution to relative power, 

as predicted by Tatarski's theoretical model. A second 
result is that for the same aperture size, the power spectrum 
relative to the power at 20 Hz measured on a day when wind 
speed was estimated to be 16 knots over the path, was overall 
higher in magnitude at all frequencies above 20 Hz than for 


the day with wind speed estimated at one knot. 


ei. Cv VS. Cur 


As stated previously, one of the goals of the project 
is to be able to predict, through meteorological measure- 
ments, optical propagation conditions. A common indicator 
of this is the refractive index structure constant, Cy : 


which is related to the log-amplitude variance through 


oy 





equations (12) and (13). Cy can be computed from the 


temperature structure constant, Co [equation (6)], through 


the relationship 


2 


= -6 P -2 2 


T 3 


where C 1s the refractive index structure constant pre- 


NT 
dicted from temperature fluctuations. Agreement between 
memstructure constant Cur obtained from long-path measure- 


ments of the log-intensity Variance with C may be taken 


NT 
as an indication that the optical fluctuations are not 
Saturated and that the point-measured temperature fluctua- 


tions are characteristic of the path. 
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Itil. EXPERIMENTAL PROCEDURE 


A. GENERAL EQUIPMENT DESCRIPTION 
i, Detectors 
To compensate for the time dependence of the 
variance, in measurement of effects of aperture Size, a 
United Detector Technology PIN-SDP photodetector (Noise 


Equivalent Power 5 x Qe" 


watts) with fixed aperture siZe 
was used as a reference. This was operated without bias 
voltage with the circuit shown in Figure 5. 

For the variable aperture detector, a TIXL69 
Silicon avalanche photodiode (operated at its optimum S/N 
reverse bias voltage of 159V) was used to provide the high 
Sensitivity at low S/N needed for measurement at very small 
aperture sizes. The characteristics of the TIXL69 are 
given in Appendix B. 

For the first series of experiments conducted in the 
enclosed corridor laboratory, both detectors operated without 
field lens systems. The field lens images the entrance 
aperture on the detector so that the light entering the 
aperture is spread over the surface of the detector. This 
avoids fluctuations due to the image wandering over the 
detector active area which may vary in sensitivity [Ref. 28]. 
For the field experiments, the following field lens systems 


were installed: 


ee, 





(a) Avalanche Detector (Figure 6). A 50 cm focal 
length objective lens with a diameter of 4.9 cm focused the 
beam onto a 1.5 cm focal length, 9.5 mm microscope objective 
lens. Between the lenses a clear glass beam splitter 
reflected a portion of the image into an eyepiece for align- 
ment purposes. The field lens magnification was designed 
so that the image of the aperture filled the detector active 


E cm*). ive og Chl. 90) A filter centered at 


meee (1.8 x 10° 
6328 A with a 77 percent peak transmission coefficient 
preceded the lens system, and an 18 leaf iris, adjustable 
from 1.5 mm to 42 mm aperture diameter, was placed in front 
of the filter. 

(b) PIN-SDP Photodiode Detector (Figure 7). This 
field lens system is fundamentally the same as that for the 
avalanche detector with the exception that a 33.3 cm focal 
length, 2.54 cm diameter objective lens and 1.8 cm focal 
length, 9 mm diameter field lens were used to cover the 
detector active area of .05 cm’, 

2. Fan-Scan Mode of Laser Operation 

A CW Radiation 5 mWatt 6328 A laser operating in the 
TEMS mode was used for the aperture averaging experiments. 
The NPS "Fan-Scan' transmitter optics developed for ship- 
to-shore measurements was used. This system broadens the 
circular beam cross-section horizontally by reflection from 
a variable-radius cylindrical mirror, then scans this ''fan" 


vertically by reflection from an oscillating output mirror. 


Thus the narrow uniform beam scans vertically twice each 
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cycle over a detector within a square field of view of order 
20 meters side at 1500 meters (variable). Ship motion and 
background compensation are provided by measuring the dif- 
ference between maximum and minimum irradiance values. 

3. Electronic Components 

Figure 8 is a block diagram of the electronic 
components used in all aperture averaging and frequency 
analysis experiments. The voltage signal from each detector 
package was sent to a PAR 113 (Princeton Applied Research) 
low noise differential pre-amplifier with maximum gain of 
10 K. The outputs from the two PAR's were connected to a 
Single switch box, from which the selected output passed 
to the demodulator. This switch box enabled the operator 
to rapidly record alternate scintillation signals from the 
two detectors. The demodulator stretched the voltage peak 
and sampled the stretched value, sampled the noise voltage 
100 microseconds later in between peaks, and sent the dif- 
ference to the Log Converter. The log-converted difference 
was then sorted by the PIP into the appropriate voltage 
channel as one count. The demodulator provides a trigger 
pulse for a signal monitoring oscilloscope and the PIP. 

The calibration box provides calibration voltages 
of 5, 1, 2, 4 and 8 volts to the Log Voltmeter/Converter, 
and from there into the appropriate voltage channel on the 
PIP 400 Multichannel Pulse Height Analyzer (MPHA). A 
readout of the counts in each channel during calibration 


is used by the HP9810A Calculator and Plotter to scale the 
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voltage. The PIP 400 MPHA records the number of times a 
particular voltage value (1i.e., log-intensity) is measured 
from the fluctuating signal during a run. This information 
is displayed during collection and can be stored indefinitely 
for subsequent printout and analysis as long as no other 
data is recorded further on that channel. The scintillation 
Signal for frequency analysis was recorded from the output 
of the Log Voltmeter/Converter. 

4. Shipboard Detector Stabilization 

The avalanche and PIN photodiode detectors were 

mounted on a feedback-controlled gyrostabilized platform. 
The system is designed to compensate for ship motion and 
keep the critical detector alignment stable. Currently this 
system's performance is adequate at short ranges (less than 
2000 m) and fairly calm sea states, but is degraded rapidly 
by increased distances or moderate sea states. In these 
instances the detector (gyroscope) alignment must be con- 
trolled manually by keeping the laser signal image in the 
center of two crosshairs of an eyepiece mounted on the 
gyro support. Gyroscope nutation, inability to keep the 
detectors aligned properly, and movements of the ship all 


contribute errors to a measured value of signal variance. 


ieee APERTURE AVERAGING EXPERIMENTS 
1. Propagation Paths and Turbulence Conditions Present 
a. Corridor Laboratory 
Aperture averaging experiments were conducted 


in a long basement corridor by locating a 5 milliwatt He-Ne 
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(6328 A) laser at one end of the corridor, propagating the 
beam the length of it, reflecting it back down the length 
of the corridor using a plane mirror and measuring the 
Seimtiliation near the location of the source. Total path 
length was 280 m. The avalanche and PIN photodiode 
detectors were clamped to a vertical Support. Once aligned, 
the detectors were not subject to any of the inherent motions 
from platform instability which are present when conducting 
the same experiment on board a ship. The experiments were 
conducted during early morning hours (with the exception 
of that on 7 April) to eliminate rapid variances in turbu- 
lence conditions due to persons opening and closing exterior 
doors. Turbulence was provided by ventilation system openings 
and building heating ducts located near the propagation path. 
Earlier measurements on the nature of the turbulence in 
this corridor indicated that the turbulence present did not 
represent homogeneous (Kolmogorov) turbulence conditions. 
However, the non-homogeneous turbulence did produce an 
intensity distribution that proved to be log-normal as 
predicted from theory. A sample log-intensity distribution 
taken during a corridor experiment is shown in Figure-9. 
b. Over Land 

To complete an aperture averaging sequence in 
homogeneous (Kolmogorov) turbulence before going to Sea, 
Where the ship's motion can affect the measured signal 
variance, an aperture averaging sequence was conducted on 


the Naval Postgraduate School Annex grounds on which the 
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Fleet Numerical Weather Central is located. The site 
selected to conduct the experiment consisted of a level 
terrain, three fourths of which was a blacktop covered road, 
the remainder a level dirt road. There were buildings and 
trees within 50 m of the propagation path; however, no 
forced draft blowers, heating systems or machinery was 
present to disturb the apparent homogeneity of the path. 
Winds were negligible; hence no eddy effects from winds 
curling around buildings or trees were present. The 5 mW 
He-Ne (6328 A) source was mounted on a large telescope 
frame in a small, open-end trailer. It transmitted a 
Fan-Scan signal over a path length of 336 m to the avalanche 
and PIN photodiode detectors mounted on a metal plate 
affixed to a manually trained telescope. The equipment 
setup of Figure 8 was used and a frequency analysis experi- 
ment was conducted in conjunction with the aperture averaging. 
Again, intensity distributions over the path proved to be 
log-normal. 
c. Over Water (Land to Land) 

Two experiments were conducted over different 
path lengths and under different turbulence conditions. 

Experiment I: Coast Guard Wharf to Spanagel 
Roof. The optical path, of length 2335 m, was highly 
irregular. The propagation path was slanted downwards, 
going from the roof at an altitude of about 75 m above sea 
level down to the Mobile Optical Research Laboratory at the 


Coast Guard Station parking lot at about 10 m above sea level. 
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About one third of the path was over water, the remainder 
over several commercial and residential structures. During 
the course of the experiment there were intermittent rain 
squalls. Data revealed that this entire Sequence was 
conducted under saturation conditions for Tor 

Experiment II: Monterey Holiday Inn to Hopkins 
Marine Station. The second experiment, conducted over a 
propagation path length of 4230 m, was 97 percent through 
the marine boundary layer. The laser beam was transmitted 
over a horizontal path at about 10 m above sea level. The 
sky was clear on 11 March with average wind speeds during 
the experiment of about 7 m/sec. The He-Ne (6328 A) laser 
source was located in the rear of the open-end trailer and 
was Operated in the Fan-Scan mode. The trailer was located 
at the Monterey Holiday Inn parking lot, about 100 m from 
the water's edge. The terrain between the ocean and the 
trailer consisted of a low, sandy incline sloping downwards 
from the trailer's position towards the ocean. The detection 
devices and processing system were located in the bus (Mobile 
Optical Research Laboratory) stationed on the Hopkins Marine 
Station grounds about 10 m above the ocean surface and 20 m 
from the water's edge with an unobstructed view of the 
Signal source. 

Experiments measuring MTF being conducted 
Simultaneously with aperture averaging precluded the PIN 
photodiode detector's being available for use as the fixed 


aperture system on this occasion. For this experiment a 
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HgCdTe detector with an aperture size of 23.6 mm was used 
as the fixed aperture system. The HgCdTe detector measured 


Scintillation from a (10.6 um) CO, laser. It is assumed 


Z 
i measured at 10.6 wm varied temporally in the same 


that C 

way as measured at .6328 um, as indicated by the consistency 

of multiwavelength measurements of ae 
d. Over Water (Land to Ship) 

The He-Ne Fan-Scan transmitter on the bus was 
located at Pt. Pinos (Figure 10) as near to the shoreline 
as possible (about 50 m away), with the detectors and data 
processing system on board the Research Vessel ACANIA. Only - 
those runs were recorded in which detector alignment was 
maintained for periods of at least, two to three minutes. 
Path lengths used during these runs were 1442 m, 1550 mn, 
and 1625 m; the maximum range was limited by the failure 
of the tracking system to lock-on automatically at ranges 
greater than 2000 m. Path lengths were determined by base- 
line triangulation methods on shore using surveyor's 
transits. The optical paths were about 98 percent over the 
ocean surface. 

ob Data Reduction 
Data reduction, i.e., determining the refractive 
index structure constant and the log-intensity variance 
from a fluctuating voltage signal, was accomplished within 
minutes of the completion of a data collection run through 


the use of the Victoreen PIP 400 Multichannel Pulse Height 


Analyzer and the Hewlett Packard 9810A Calculator and 
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Plotter. As described previously, the PIP 400 stored the 
log-intensity distribution of the irradiance measured. 
An interface was developed between the PIP 400 and the 
HP9810A to allow the calculator/plotter to print out the 
from 


log-intensity distribution and determine C,, and o 


N RI 
it. The algorithm the computer followed [Ref. 31] was to 
apply a best-fit Gaussian curve to the distribution to 
through the relation 


obtain o and then determine C 


XT? N 


Mi 


‘il eee 2 Tz 
Nk Z (.476) 


An example of the printout is shown in Figure 9. The three 
markers near channels 60, 116 and 173 are calibration 
points of .5 V, 2 V and 8 V respectively, which are equally 
spaced on the log scale. 
3. Aperture Averaging Factor G(D) 
a. Description 
To determine the aperture Size that best 

approximates a point receiver, l.e., that size at which 


Oo reaches a maximum for a given path length and Cx 


QI 
value and does not decrease significantly for further 
reduction of aperture size, the following aperture averaging 


factor was used: 
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Cy[V(D,.t,)] CylF(D).t)] Cy(V@,,t)] 
(48)  G(D) = . : : 
mi G@emcaie CIV, cj]  C.1vOD,t)] 


Cy [V(D,>t,)] is the value of Cy measured using the variable 
aperture detector at the maximum diameter De and the 

time ty: Similarly, F denotes the fixed aperture detector 
which is used to measure the change in Cy due to variables 
meter than a change in aperture size. The first factor 

Cy fV(D,>t,) 1/Cy [F(D,, tJ], compares the initial Cy value 
measured (at time cae with the variable diameter set at . 
to that measured by the fixed diameter at time t.-: 

Cy [F(D,,t) 1/C, (VD, t)] reflects the change in that ratio 
produced by a reduction in aperture diameter from Dy COr2 


If the two apertures (fixed and variable) measure o Simul - 


1 
taneously at times to and t, and the variable aperture size 
is not changed between the two measurements, then G(D) 
Will ideally equal 1.00. This was tested in conjunction 
With two corridor experiments conducted on 10 February 
and 5 March 1976. The comparison runs produced G(D) values 
of .9818 and .9981, indicating an average error of about 
one percent for G(D) values at different aperture sizes 
due to experimental variances. 
b. Reason for Use; Relation to Fried's 0O(D) 

By starting with the largest experimental 

variable aperture diameter, Da» for which G(D,) = ecie. 


then progressively reducing it until the minimum variable 


48 





aperture Di is reached, the aperture size that best approxi- 
mates a point detector can be determined by noting at which 

aperture diameter the lowest value of G(D) occurs. Empirical 
results (Section IV.A.) show that G(D) does not remain linear 
from D, down to the smallest possible aperture size, Din? 


but levels off and reaches a minimum value, [G(D)] well 


mene 
before it. 

Fried [Ref. 21] defined an aperture averaging 
factor O(D) [equivalent to A in equation (41)] based upon 
a vanishingly small receiver (D + 0 mm). Fried's theory 
States that all apertures of diameter D > 0 will measure a 
Signal variance of smaller magnitude than would be measured 
using a vanishingly small detector. When using an aperture 
of diameter D, the measured variance must be multiplied by 
the correction factor — to allow for aperture averaging. 


Rod) in identifies the aperture diameter, D , that best 


min 
approximates a vanishingly small receiver. G(D) is related 


to Fried's aperture averaging factor in the following manner: 


4g) 8D) = CNEVOg2 8] CyIVD,8)] _ CylVOQ4] 
Ran Cy IVT, £7] yl oe ey VOR n? 2] 9(D) 


C. APERTURE AVERAGING EFFECTS ON TEMPORAL-FREQUENCY 
SPECTRUM ANALYSIS 


1. Experimental Conditions 
a. (7-9) April 1976, Basement Corridor 
The electronic equipment setup used in aperture 


averaging sequences is shown in Figure 8 and the procedure 
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followed in obtaining and processing signals on this 
date is described in Sections III.B.l.a, IV.A.1, and 
mex. 2.a.(3). For frequency analysis, the scintillation 
Signal was recorded from the output of the Log Voltmeter/ 
Converter with a Precision Instruments, Model No. PI-6204, 
four-channel recorder at 37.5 ips. Performance characteris- 
tics of the recorder are contained in Appendix C. Signals 
were recorded for a period of about 75 seconds each at 
eeemecune Sizes of 40, 30, 20, 15, 10, 7, 5, 3, 2 and 1.5 
millimeters. 
b. 20 April 1976, Navy Annex 

Experimental conditions are as described in 
Beectons L[If.B.1.b, IV.A.1, and IV.A.2.c. Frequency analysis 
data was recorded as described above for the (7-9) eel EUs. 
Aperture diameters examined were of 40, 30, 20, 15, 10, 7, 
5 and 3 millimeters. Power to operate the equipment on 
the bus was provided by a portable gasoline-powered motor 
generator. 

c. (29-30) April 1976, R/V ACANIA 

Experimental conditions on this date and 
procedures followed for aperture averaging are described 
mamoeeti1ons I11,B.1.d, IV.A.1, and 1V.A.3.(d and e). The 
Signal to be frequency analyzed was recorded from the 
output of the Log Voltmeter/Converter (Figure 8) on a 
Sanborn HP-3914, 14-channel Analog Tape Recorder, operated 


by the meteorology team recording micrometeorological data 
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for Cy comparison. Operating characteristics of the 
HP-3914 Recorder are included in Appendix C. Time reference 
Signals were recorded on the tape for identification. 

2. Data Reduction 

The recorded scintillation signal was frequency 
analyzed by playing back the signal into an EMR Schlumberger 
Instruments 1510 Digital Spectrum Analyzer. The spectrum 
analyzer was operated in the continuous mode over a frequency 
range of 10 Hz to 2.56 KHz. The signal was sampled by the 
Spectrum analyzer for a period of 25.6 seconds. The resulting 
Spectrum was plotted on a Hewlett Packard Moseley 7035B X-Y 
Recorder. Frequency was plotted linearly along the X-axis 
from 42 Hz to 2.56 KHz. The Y-axis was a 0 to -60 dB linear 
sicaie, 

For evaluation of the aperture averaging effects on 
the spectral density, relative power spectra were obtained 
directly by aligning the -60 dB line of the spectrum of each 
different aperture to a common new reference axis and 
plotting each spectrum manually. Points of local minima 
were plotted on each spectrum, thereby avoiding abnormal 


magnitudes due to noise spikes. 
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IV. EXPERIMENTAL RESULTS 


A. APERTURE AVERAGING EFFECTS ON SCINTILLATION 
1. Commonalities 
In addition to the electronic equipment setup 
(Figure 8), detectors, and laser source previously described, 
the following items were common to all aperture averaging 
experiments: (1) The Hewlett Packard 9810A Calculator and 
Plotter was used for immediate data reduction. (2) Aperture 
Sizes for the variable aperture system were set by inserting 
a micrometer caliper set at the proper diameter size within 
the aperture and closing the iris upon it. Both horizontal 
and vertical dimensions were checked. The aperture size was 
rechecked at the end of the run. 
2. Preliminary Investigations 
a. Corridor Experiments 
(1) 2235, 7 February 1976 - 0515, 8 February 1976 
This initial aperture averaging sequence was 

conducted to examine the validity of the assertions made 
by Schroeder [Ref. 28] that a 4.8 mm aperture diameter best 
approximates a vanishingly small receiver. For this experi- 
ment, the laser was operated in the Scan mode only and a 
neutral density filter of .5 was used to decrease the signal 
Strength. The field lens detector systems of Figures 6 and 
7 had not been constructed at this time. Each detector 


consisted of the photodiode placed at the focal plane of an 
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objective lens in a narrow cylinder of diameter 22.6 mm with 
a 10 A filter centered at 6328 A at the cylinder front. This 
made the largest variable aperture size available 22.6 mm. 
M@erture sizes of 22.6, 20, 17, 14, 12, 10, 7, 6, 5, 4, 3 

and 2 millimeters were used. G(D) values calculated from 

the measured log-intensity distribution standard deviations, 
Soy, are given in Table I with the plotted "G(D) vs. Aperture 
Diameter, D'" shown in Figure 11. The last column of Table I 
G(D) 
D 


, the correction factor for aperture 
min 


gives ann 
averaging. 
In Figure 11, the curve bottoms out at an 
aperture size of about 5 mm and the corresponding value of 
G(D) does not significantly decrease from its value at 
that point as the aperture diameter is further reduced. 
(2) 1655-1820 and 2100-2240, 10 February 1976 
It was noted from the experiment of 
7 February that although the log-intensity distributions 
were for the most part Gaussian, on several of the larger 
aperture runs double peaks occurred near the top of the 
Peetncensity curve. lt was believed that this was due to 
the laser being in the Scan mode only so that the detector 
would see one average intensity peak as the beam swept up 
and a slightly different one as the beam swept down. In an 
attempt to alleviate this, a second aperture averaging 
Sequence was conducted with the laser in the Fan-Scan mode. 
By spreading the beam out and scanning it, it was felt that 


a more uniform intensity might be seen by the detector, and 
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the double peaks eliminated. <A Fan-Scan divergence setting 
of seven milliradians was used and a neutral density filter 
of optical density .5 was inserted in the beam to reduce the 
Signal intensity at the detector. With the beam reflected 
Set a 625 cm? mirror at a distance of 140 m, the estimated 
Signal strength at the position of the detector is 1.7 x 107° 
Watts/em’, which is of the order of the signal intensity 

for the at-sea experiments, described later, at a path 

length of 1600 mn. 

Aperture sizes used were identical to those 
used during the 7 February run. Log Intensity distributions 
on this day revealed that the double peaks were eliminated 
in all but one distribution by use of the Fan-Scan, and 
this mode was therefore adopted as the standard laser opera- 
tion mode for future aperture averaging sequences. Figure 
12 shows "G(D) vs. Aperture Diameter, D'" for this experiment. 
Although this experiment was conducted identically to that 
on 7 February with the exception of the Fan-Scan operation, 
three significant differences stand out in the results. 
First, the curve bottoms out at the 6 mm aperture size vice 
the 5S mm size of 7 February. Secondly, G(D) values at all 
aperture sizes are in general larger than they were on 
7 February, even though the same value of D, (22.6 mm) was 
used to calculate G(D) values. Thirdly, the shape of the 
curve following the bottoming out at 6 mm is different from 


that seen on 7 February. In comparing the data, the 
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average o,,[F(D,>t)] Hemceen to be Signiticantly different 


on the two days, implying a larger C On 7 February 


N° 
<do7[F(D,,t) ]> = .330, while on 10 February <o, ,[F(D,,t) ]> =.428. 
An increasing <O,7? Seems to increase G(D) values, as can be 
seen by comparing Table II with Table I. 
(3) (7-9) April 1976 

An aperture averaging frequency analysis 
experiment was conducted. The aperture averaging results 
are included for comparison with runs conducted on 7 February 
and 10 February. The same He-Ne (6328 A) laser was 
operated in the Fan-Scan mode with a divergence setting of 
seven milliradians and a neutral density filter of .5 was 
inserted in the beam. 

Prior to this run the largest aperture size 
available on the variable aperture detector was 22.6 mm. 
To increase the largest aperture size available and improve 
Signal detection capability, the field lens systems of 
Figures 6 and 7 had been constructed and were in use for 
this experiment. The significance here is that with a 
larger Da available, one would expect to see a smaller 
BED) } an for this sequence. This proved accurate, as 


Shown in Figure 13. A much lower [G(D) ] resulted using 


min 
the 40 mm aperture as Do Vice 22.6 mm. This curve appeared 
to follow the same basic pattern as those for 7 February 
and 10 February with the exception that this took on a 


different shape from either of the previous two for apertures 


less than the bottoming out diameter of about 5 mm. In 
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order to see better how this run compared with the previous 
two, G(D) values were recomputed using the 20 mm diameter 

as Do. fice rcciwmr ts Shown in Figure 14. [G(D) Js. for 
(7-9) April is still well below that of 7 and 10 February. 
The average Tor for the fixed aperture system for this run 
was <o),[F(D,,t) 1? = .08 which is much lower than the average 


values for 7 and 10 February as is [G(D)] Table III 


min” 
gives a comparison of the average signal variance with 

ree an for the three experiments conducted in the same 
manner, with the same equipment, but with different turbu- 


lence conditions present as indicated by <o,,> 





TABLES TE 
Date <o, [F(D,,t)]> ae oa Oe 
10 Feb ~428 me -- -- 
7 Feb ~ 3350 - 468 .098 .179 
(7-9) Apr .084 w2Lo ~ 246 ~245 


It appears that [G(D) ] increases with <O 97? for low values 


min 
‘aye <Ooy? > and that in general, for a given D> the value of 
[G(D) ] 


higher the turbulence level, the higher the [G(D)] 


ee is dependent upon the level of turbulence. The 


se value. 
b. Aperture Diameter Normalization 

Examination of Figure 14 shows that the shapes 

and slopes of the plots are similar from the largest aper- 


ture size down to about 5 mm where three different forms of 
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the curve then emerge. Following previous theoretical 


models [Refs. 1, 20 and 21], it is expected that D iene 


min’ 
aperture best approximating the "vanishingly small" receiver, 
should be nec aley related to the intensity transverse 
correlation length, Po. This quantity is given in the 

three models as: 


4) a2z)1/* | pratarski, Ref. 20] 


3 


(42) 1/2 


(2) (= , [Fried, Ref. 21] 


(3) (E) 1/2 Methcominns <lmand titra.) Ret gel). 
Aperture diameters for the three corridor experiments were 
scaled according to the three models, the results of which 
are shown in Figures 15, 16 and 17, respectively. The 
three curves bottom out at approximately the value of 


1/2 
D = (5) predicted by Lutomirski and Yura. 


c. Over Land Experiment (2128-2347, 20 April 1976) 

Aperture diameters examined were 40, 30, 20, 
15, 10, 7, 5 and 3 millimeters. The Fan-Scan was set at a 
divergence of 7 milliradians and no neutral density filter 
was used to decrease the signal intensity. The laser was 
situated about 2 meters and the detectors 2.5 meters off 
the ground. The sky was clear, temperature averaged 10°C 
and winds were negligible (0-2 m/sec). 

me CD) vst ere WasmeeloOgEed to Locace the 
bottoming out point. Although able to test only a few 


aperture sizes, previous runs conducted outside the corridor 
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laboratory (discussed in section d below) indicated that the 
72 
(z/k) 


the bottoming out point of the graph at D/(Z/k) 


normalization factor was accurate in predicting 


V/2 2 100. 


Figure 18 shows the graph for the experiment conducted on 


20 April. The bottoming out occurs for the aperture size 


nearest, but not less than, D/ (z/k) 1/2 


le 


= 1.0. Note, for 


this path length (Z/k) 5.8 mm, and although the 5 mm 


aperture size is nearer to (z/k) 1/2 


than the 7 mm size, 
previous runs have shown that the behavior of G(D) below 


may)’ 2 


to the minimum G(D) value that corresponds to the bottoming 
WZ 


= 1.0i1s erratic. Thus, the 7 mm size is nearer 
ome point at D = (Z/k) than is the 5 mm size. Asa 

Mercier Comparison, the corridor runs conducted on 

7 February and (7-9) April at a path length of 280 m were 
plomered With this nun of 20 April at the path length of 336 m, 
shown in Figure 19. It appears that the bottoming out point 


oe eal.) 


occurs very close to D/(Z/k) 
[G(D) J in for this experiment, with respect to 


a 20 mm D,, was .666 and the average o was 


i 
<o, 7 [F(D,>t) 1? = ,162. This does not fit in with the results 
maven in fable Ill; however, the propagation path and turbu- 


lence conditions were different. To verify a predictable 


Wath <0 42> would require repetition 


increase of rei haee aI 


of the experiment under identical experimental conditions 
and over the same propagation path where the turbulence 
would have the same characteristics. Time considerations 


precluded this procedure. 
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d. Over Water (Land to Land) 


CO 22 0eeselarch = 0345.56 March 1976 (Spanagel 
Roof to Coast Guard Parking Lot) 


Due to the. irregular propagation path used 
and corresponding turbulence conditions over a long distance, 
this entire sequence was conducted under saturation con- 
@mcions for Sor: so, ,(F@.,t)]> was 1.01 and [S(D) Jn 
was .8808 with respect to a Do Gre2¢4.5 mm, The graph of 
"G(D) vs. Aperture Diameter, D'' is shown in Figure 20. 


For this path length, (Z/k) 1/2 


= 15.3 mm. The main item 

of interest for this run is that, under saturation conditions, 
the maximum difference in G(D) values observed was only 

et D) ax: .167, which included comparisons of aperture 
peegeseo: 24.5, 20, 16, 12, 7, 6, 5, 4, 3 and 2 millimeters. 
This indicates that aperture averaging has only a small 

effect on the measured signal variance under saturation 


conditions for aperture sizes less than 24.3 mm. 


(2) 1550-1640, 11 March 1976 (Monterey Holiday 
Inn to Hopkins Marine Laboratory) 


Due to the limited time available for 
eenducting the experiment On this day, only aperture sizes 
of 41.7, 24.3 and 14.0 millimeters were examined. A plot 
of "G(D) vs. Aperture Diameter, D'' is shown in Figure 21. 


EZ ye 20.6 mm, 


The path length was 4230 m, giving a (Z/k) 
Figure 21 contains insufficient data to draw definite 
conclusions. The same data is plotted again in Figure 24 


in terms of the normalized diameter, and will be discussed 
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at that point. [G(D) ] sn and <o, [F(D,,t)]> Om eeneese Ui 
should not be compared with other experiments in this sec- 
tion since a different source and fixed aperture detector 
system was used Bed Tor 1s wavelength dependent for identical 
meppuience conditions. 


3. Propagation in the Marine Boundary Layer (R/V ACANIA 
Experiments) 


a. 1840-1930, 30 March 1976 

The winds during the time of this experiment 
varied from 8 to 12 m/sec in a direction of about 10 degrees 
off the normal to the propagation path. The sky was partly 
cloudy, temperature about 7°C, sea conditions were moderate 
With average wave heights of about 1 meter, and swell heights 
of about 2 to 3 meters moving in a direction normal to the 
propagation path. The motion of the R/V ACANIA at anchor 
made maintenance of gyro lock-on difficult. This undoubtedly 
caused intensity variances to be greater than would be ob- 
served normally. The path length on this day was 1600 m. 
The detectors were located about 7 meters and the source 
about 5 meters above the ocean surface. Aperture sizes 
ecompared were 38, 31, 24.3, 17 and 10 millimeters. For 


mis Pat length, (z/k) 1/2 


zie mm, Lhe e(D) vs. 
Aperture Diameter, D" curve is shown in Figure 22. This 
1s not typical of past aperture averaging results, but 
resembles those obtained under saturated conditions. How- 


ever, <o 7 LF(D,,t)]> for this experiment was .149, well 


below the saturation region. It is felt that the sea 
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Somditions , rapid ship movements, and inability of the gyro- 
scope lock-on system to work effectively in moderate seas 
may have obscured any detectable aperture averaging effects. 
b. 1340-1415, 2 April 1976 

Wind and sea conditions on this day were 
identical to those of 30 March. The sky being clear was 
the only difference in weather conditions. The sea state 
was Still moderate and there was still considerable ship 
motion. The path length on this day Teeeelenn m. The gyro- 
scope lock-on system had been improved slightly by readjusting 
the existing circuitry; all other experimental conditions 
were as described for 30 March. Aperture sizes used for 
this run were 40, 30, 20 and 10 mm. Figure 23 shows 
Zoe vs. Aperture Diameter, D" and Figure 24 compares 


1/2 


mo) vs. D0/(Z/k) " for this run with the data taken over 


water with a path length of 4229 m on 11 March. Once again, 


1/2 successfully 


it appears the normalization factor (Z/k) 
predicts the diameter of the aperture that best approximates 
a vanishingly small receiver. 
Ceeelad o—ooy 27 April 1976 

The sky was clear at 1445 but became overcast 
at about 1545. Winds were 8 to 10 m/sec blowing normal 
to the propagation path; waves also moved normal to the path 
With a height of about 1 meter and swells about 2 meters 
high. The ship rolled considerably in the swells and once 


again gyroscope lock-on was difficult. The ship was 


anchored, path length was 1650 m, and aperture diameters 
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memmnecd were 405655, 50, 25, 20, 16, 13, 10, 5 and 4 


teplameters, "G{(D) vs. D/ (2/k)t/? 


" for this run is shown in 
Pagure 25. G(D) values seem erratic; this is attributed to 
difficulty in stabilizing the detector and tracking the 
ACANIA from shore. 
d. 1830-2024, 29 April 1976 

The sky was clear, wind speed throughout the 
experiment was 3 to 4 m/sec, and path length was reduced 
to 1442 m. Sea conditions improved slightly with wave 
heights down to about .7 m and swells of about 1.5 m. The 
wave direction was about 15 degrees off the normal to the 
propagation path as was the wind. Aperture sizes of 40, 35 
30, 25, 20, 16 and 13 millimeters were used. "G(D) vs. 
D/ (2/k) 14? 


Size was used as Do for the G(D) calculations since it 


" is plotted in Figure 26. The 35 mm aperture 


appeared to give more credible overall values based on 
past experiments. The G(D) value for the 40 mm aperture 


1s plotted relative to the 35 mm aperture as Do: With the 
D 


exception of the 40 mm aperture eee = 3.33] and the 
(Z/k) 
D 
16 mm aperture ie 1.33], the curve appears to 
(Z/k) 


parallel results of past runs, i.e., reaches a minimum at 


7” 1.0 and shows a nearly linear decrease of G(D) 
(Z/k) 


With decreasing aperture diameter to this point. For this 


mim, [G@eD) } = .721 and <o,,[F(D,,t)]> = .145. 


min 
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cela ollie oO April 1976 

The sea state was relatively calm with wave 
heights of about .5 m and no apparent swells. The sky was 
elear, Meera cure about 21°C and winds were 3 to 4 m/sec. 
The calm sea state enabled the gyroscope tracking system to 
work well in the automatic mode and results more consistent 
with past sequences were obtained on this day. The ship was 
at anchor and the propagation path length was 1625 m. 
Aperture sizes examined were 40, 16, 13, 10 and 5 milli- 


meters and (z/k) 1/2 


for this path length was 12.8 mn. 
The plot is shown in Figure 27. Bottoming out occurs near 


D/ (Z/k) 1/2 


<o,7[F(D,,t)]? = ,264. <A comparison of [G(D)] 


= mela [G(D) J win =eenoo wand Lor thiserun 


- and 
min 


<O 97? for the three R/V ACANIA experiments is given in 





Table IV. 
TABLE IV 
Date [G(D) ] an <o,,[F(D,,t)]? 
2 Apr De eco 
30 Apr mie . 264 
29 Apr eet ~ 145 


From the data of Table IV it appears again that in general 


[G(D)] 


; increases with <o,-+> 
ae ncrease t at 
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4. Comparison of Theoretical Aperture Averaging Factors 
Wit Mpc tiene a SSulesS 


a. Lutomirski and Yura Model [Ref. 1] 
For aperture diameters on the order of the 
intensity correlation length, as utilized during these 
experiments, Lutomirski and Yura "recommend" an aperture 


averaging factor of 


(46) ys 


Table V lists Lutomirski and Yura's theoretical value of A 
computed for our experiments and the observed value of A, 


A, theoretically defined as s/o. ‘ , 1S equivalent to 
O 


our [G(D)] /[G(D,)] where D, replaces D in equation (46). 


min 
The theoretical values do not correlate well with the 
empirical values. Further, Lutomirski and Yura's theoretical 


model does not contain any turbulence dependence which appears 


to be necessary in view of the apparent dependence of 


[GCD) ] Seo ccen cmp iuteally.: 


min QI 
b. Fried's Model [Ref. 21] 


To compare empirical results with Fried's 


theoretically derived plot of "0 vs. D/ (4Z/k) 1/2" 


(Figure 4), whose curves are based on an aperture diameter 


2 


oe = 2-1(42/k} as representing a vanishingly small 


receiver, it was necessary to use the experimental aperture 


66 





LoS 
Og Si 
Loo. 
v8c° 
979° 
ECOG 


612° 


CIVOIUIdND) V 


062° 


Coe 


882° 


LV oe 


S6v- 


990° 


OSt: 


(CRIOFHL) V 


A HTave 


VINVOV/idy 0¢ 
VINVOV/1dy 62 
VINVOV/idy Zz 

xouuy AAen/ady QZ 

APT TOH- sutydoy/i1eW TT 
Loptiioj/ady (6-Z) 


LOpTiI10)/qeyq 1 


NOILVIO1T/SLVG 


67 





Size for each run that gave the minimum G(D) value as the 
vanishingly small aperture diameter. The effect of this is 
to shift the empirical curves to the right of the theoretical 
curve. Figures 28 and 29 show Fried's theoretical plot 

for log amplitude variance C, (0) = ae as compared with 
experimental values for various runs. For this graph, 
aperture diameters were appropriately normalized for com- 


ae 


MPaeison by the (4Z/k) factor. Using equation (45), 


0 was calculated as 


Z 
o."[V(D,t,)} 


(47) @ = =, 
omens) 


where o.“[V(D',t.)] is the etearee variance measured at 

time Ee and with variable aperture diameter D', which 
corresponds to the aperture diameter that best approximated 
a vanishingly small receiver, i.e., corresponded to 

[G(D) ] win: o.*[V(D,t)] 1s the signal variance that would 
have been measured at time ey had aperture size D been used, 


and was found from the relation 


Z 
cei Oe ie 


Z 2 
(48) 9. “(V(D,t,)] = 9° (V@,t)] - 4, 
: ° = o.[F(D,>t)] 


where again, F is the fixed aperture detector. 
Empirical results appear to follow Fried's 


theoretical model generally well. In contrast with 
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Lutomirski and Yura's aperture averaging model, Fried 
acknowledges the dependence of the aperture averaging factor 
on the turbulence conditions, as represented by different 
curves for pereuienilar values of C, (0) Pie Gourertae 1hied's 
model indicates that as C, (0) increases, 98 increases for a 
Specific aperture diameter, or the inverse of the aperture 
averaging factor decreases. This 1S in agreement with 
empirical results noted in Sections 2.a.(3) and 3.e. above. 
From Figure 28 it is evident that the aperture averaging 
curve for 7 February lies below that of 10 February and above 
that of (7-9) April. The 98 values are therefore in 
general larger for 10 February than for 7 February, il.e., 


the ratio of o.*/s. 2 
O 


is larger, and the inverse of the 


aperture averaging factor is less, as predicted by Fried's 
theory. 
B. APERTURE AVERAGING EFFECTS ON THE RELATIVE TEMPORAL- 

FREQUENCY POWER SPECTRUM 

1. General Results 

The spectrum analyzer frequency limits were selected 

to examine frequencies within the range 10 Hz to 2.56 KHz. 
The demodulator in the signal processing system eliminates 
frequencies above 2 KHz, coinciding with the original 
estimate of the Fan-Scan frequency. The bandwidth is further 
reduced to a maximum frequency of 1 KHz for experiments 
conducted using the Sanborn HP-3914 Recorder. Power 


spectra from all frequency analysis experiments show the 
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actual Fan-Scan frequency to be 1.92 KHz, clearly evident 
as a noise spike in Figures 34, 35 and 36. Previous work 
at NPS by Schroeder [Ref. 28] examined frequencies from 
mdetiz to 1 KHz. } 

Results of the experiments conducted on (7-9) April, 
20 April and (29-30) April, under varying turbulence 
conditions, all indicate that aperture averaging has a 
discernable effect on the observed power spectrum. Empiri- 
cal results support the theoretical prediction that a 
corresponding reduction in high frequency power results as 
aperture diameter is increased. The reduction in power, 
however, appears to occur over a large portion of the fre- 
quency range, not Solely on the high frequency end. Figures 
30, 31, 32 and 33 show the effect of aperture averaging on 
the relative power spectrum. 

2. Turbulence-Dependent Forms of the Frequency Spectrum 

As previously remarked, the three frequency analysis 
experiments were conducted at different locations and under 
different turbulence conditions. The resulting power 
spectra for each run were similarly different. Typical 
spectra for the corridor experiment (inhomogeneous turbulence), 
the Navy Annex run (over land in homogeneous turbulence), and 
the R/V ACANIA data (over water in homogeneous turbulence), 
are shown in Figures 34, 35 and 36, respectively. The 
Significant item of note in Figure 34 is the uniformity of 
the spectrum below .85 KHz. The prevalence of noise spikes 


at 110 Hz and 120 Hz intervals in Figure 35 is believed due 
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to the extremely noisy, gas-powered portable generator 
used to provide power to the electronic equipment during 
this experiment. Also peculiar to this experiment were 
@iaracteristically low magnitude power spectra. Figure 36 
Shows an approximately linear decrease in power with 
increasing frequency. 

Eliminating the noise spikes in Figure 35, it 
appears that the spectrum assumes the same general form as 
that of Figure 34. A significant difference in the experi- 
mental conditions for 9 April and 20 April as opposed to 
those of 30 April was the average wind speed present. 

On 9 April and 20 April winds were negligible. On 30 April, 
conducted at sea on board the ACANIA, average wind speeds 
were about 4 m/sec. Further experiments need to be con- 
ducted to examine the effect of wind speed on the shape of 
the power spectrum. 

3. Comparison with Clifford's Theoretical Model 

Clifford's theoretical model for the temporal- 
frequency power spectrum for spherical and plane waves 
are given by equations (32) through (35). The values of Ce 
for 9 April, 20 April and 30 April, based upon respective 
average wind speeds of 1 m/sec, 1.5 m/sec and 4 m/sec are 
30 Hz, 41 Hz and 53 Hz, respectively. For comparison with 
the theoretical model (Figure 3), typical spectra for these 
experiments were plottedvs. "Log 2" and are shown in 


Figure 37. The vertical axis is the percentage of power 
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at a particular frequency compared to the maximum power 
attainable as determined by the spectrum analyzer's 
"Y-axis" voltage setting. 

There rear to be two significant differences 
between the empirical results of Figure 37 and Clifford's 
m@ieonetical spectra of Figure 3. First, Clifford indicates 
that the rapid downward slope of the spectrum begins at 
log 2 *-0.1 for plane waves and log Q2 * 0.2 for spherical 
waves. Empirical results show the drop occurring at 
log 2 * 3.0 for the (7-9) April and 20 April experiments. 
The spectrum of (29-30) April shows an overall negative 
Slope for all frequencies. Secondly, the slope appears to 


-8/3 


decrease less rapidly than the theoretical (f) ihe 


Gy del 


Slope for (7-9) April is of the order of (f) and for 


2/13) Further data 


Homeport) and (29-30) April, (f£) 
collection is necessary to examine if, and to what extent, 
various turbulence parameters affect the magnitude of the 
negative slope. It does appear, however, that Clifford's 
(£) 8/5 1s too rapid a decrease. 


4, Differences in Magnitude of the Aperture Averaging 
EPree t 


Examination of Figures 30 through 33 indicates 
that there are significant differences between aperture 
averaging effects on the relative power spectrum as seen 
under the various experimental conditions. In an attempt 
to determine a factor common to the variances and to the 


power reductions which were not directly proportional to 


ez 





aperture diameter, parameters such as wind speed, <o,;> , 
Cy Oot and G(D) were examined. No factor was found 

that individually accounts for the variations. A detailed 
analysis of all factors concurrently was not undertaken 


due to time limitations. 
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V. SUMMARY AND CONCLUSIONS 


A. APERTURE AVERAGING EFFECTS ON SCINTILLATION 

The variance of log-amplitude fluctuations that would 
theoretically be observed if using a vanishingly small 
receiver can best be approximated by an aperture of diameter 
oe! 2. 
averaging effects present when using apertures larger than 
(z/k) 1/2, 


diameters smaller than (Z/k) 


A diameter of this size eliminates aperture 


and that may be present when using aperture 
py Experiments conducted 
in an enclosed laboratory in inhomogeneous turbulence and 
field experiments conducted in homogeneous turbulence, 
both over land and in the marine boundary layer, support 
this prediction from the theoretical model of Lutomirski 
and Yura [Ref. 1]. 

Fried's theoretical model of aperture averaging 
[Ref. 21] predicts that 9, the aperture averaging factor, 
increases continuously as aperture size is decreased from 
an arbitrarily large aperture diameter. Empirical results 
indicate that the aperture averaging factor does not increase 
continuously as aperture size is reduced to a vanishingly 
small diameter, but increases continuously until an aperture 


1/2 impevened. A L£UEther Feduction in 


eZ 


diameter (Z/k) 
aperture diameter [below (Z/k) ] causes 9 to (1) reverse 
direction and decrease or (2) remain near the same value 


Ou” . 
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In general, the value of the aperture averaging factor 
increases as <0), increases, Put another way, for an 


aperture of diameter larger than (z/k)1/? 


, the aperture 
averaging correction factor decreases as the amount of 
turbulence increases. Under weak turbulence conditions, 
the effect of aperture averaging is quite significant. For 
Strong scintillation conditions, the aperture averaging 
effect is decreased and becomes less significant. 


The ratio G(D')/[G(D)] gives the factor by which the 


min 
measured signal variance (or other turbulence parameter), 
a2 


3 


measured with an aperture of diameter D', D' > (Z/k) 
must be multiplied to represent the value that would have 
been measured had a vanishingly small receiver been used. 
It is important to note that the value of Ee an is a 
function of the largest aperture diameter used (DJ; and 
should not be viewed as an absolute value for a particular 
set of turbulence parameters. In general, the value of 


[G(D) ] 


. decreases as D_is increased. 
min O 


Pee DekwuRE AVERAGING EFFECTS ON THE TEMPORAL- FREQUENCY 
POWER SPECTRUM 
An increase in Bet rUre diameter appears to decrease 
the magnitude of power in the frequency range 
.5 KHz < £ < 2 KHz of the temporal-frequency power 
spectrum for scintillation measurements. The difference 
in power at a particular frequency for two given aperture 
diameters differed significantly for different experimental 
and turbulence conditions. The reason for the disparity 


Was not determined due to time considerations. 
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Spectra obtained from experiments in the corridor 
laboratory, over land, and in the marine boundary surface 
layer indicate a negative sloped decrease in power at 
higher frequencies for the temporal-frequency power 


spectrum less rapid than the (£) 78/5 


predicted by 
Clifford [Ref. 19]. Further investigation is required 
to determine the effect changes in certain turbulence 


parameters have on the form of the power spectrum. 


Gee eXtENSION OF PROJECT 


Comparisons of [G(D)] With <o,.> indicated that, 


we 


for a given aperture size larger than that which approxi- 


min 


mates a vanishingly small receiver, the aperture averaging 
correction factor 1s not a constant, but is dependent upon 
the level of turbulence present. Further experimentation 


is needed to compare [G(D) ] values, for a given maximum 


min 


reference diameter D> With <o,,> values, to determine the 


al 

relationship between the level of turbulence and the aperture 

averaging correction factor for a given aperture diameter D,: 
It seems clear from this report and other published 

data that aperture averaging reduces the high frequency 

power content in the temporal-frequency power spectrum 

and eliminates certain high frequencies proportional to 

increasing aperture diameter. However, experiments conducted 

under different turbulence conditions and over different 


propagation paths indicate that the power spectrum assumes 


different general shapes. This needs further investigation 
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through the examination of power spectra obtained over the 
Same propagation path under various turbulence conditions, 
and under the Same turbulence conditions Over various 
Meepagation paths. In conjunction with this, Clifford's 
theoretical model of the temporal-frequency power spectrum 
for plane and spherical waves should be examined for validity 


more thoroughly through empirical study. 
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Figure 1 


Illustration of Kolmogorov Turbulence Model 
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Figure 3 


SPHERICAL WAVE 


-0.8 ° -0.4 0 0.4 25 


log & 


Clifford's theoretical temporal power spectrum of log- 
amplitude fluctuations for plane- and spherical-wave 
propagation. 


82 





Figure 4 
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Figure 8 


Electronic Components: Aperture Averaging Experiment 
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APPENDIX B 
TEST PROCEDURE AND PERFORMANCE CHARACTERISTICS 
OF THE TEXAS INSTRUMENTS TYPE TIXL69 LARGE-AREA 
SILICON AVALANCHE PHOTODIODE 
A series of tests was conducted to evaluate the 

detection capabilities and performance characteristics of 
the TIXL69 avalanche photodiode detector and to compare 
them with that of the PIN-SDP photodiode detector. The 
two detectors were evaluated at a Sees ee rere compara- 
ble to that anticipated during a typical at-sea data run, 


i.e., about 1.4X107°W/em? 


(the intensity of a 5 mW He-Ne 
(6328A) laser source at a distance of about .8 nautical 
miles and with a beam divergence of 7 milliradians). 

The circuit used with the avalanche detector is 
illustrated in Figure Bl(a). The reverse bias voltage 
source was a Hewlett Packard Harrison 6515A DC Power 
Supply. 

Dark current as a function of reverse bias voltage 


was measured first to see if it compared with the manu- 


, 7 


facturer's stated values of 3.5X10 “A minimum to 1.1X10 ‘A 
maximum, The test circuit used is shown in Figure Bato). 
During this experiment the avalanche detector was covered 
Euebiackecloth to keep light from hitting the detector 
surface. There were no windows in the enclosed room, all 


room lights were out, and the room was completely dark. 


A flashlight was used for recording data. To ensure that 


SLL 





the voltage applied was actually occurring as a voltage 

drop across the avalanche photodiode, a Hewlett Packard 

419A Null Voltmeter was connected to the circuit as shown. 
The voltmeter imaiea ted that the avalanche photodiode was 
indeed receiving the full potential drop. Using a Keithley 
Phiesmneuments 6008 Electrometer to measure the current, it 
appeared that values in the nanoamp region were occurring. 
However, the current indicator was not stable and repetitive 
current readings at the same reverse bias voltage fluctuated 
greatly. It was decided that a more accurate ammeter was 
needed and a Keithley Instruments 410 Micro-Microammeter was 
obtained. Steady, repeatable values of current were now 
attained for specific reverse bias voltages. 

As shown in Figure B2, dark current was linear with 
applied reverse bias voltage from 20 volts through 157.5 
volts, then increased by a factor of about Loe at an 
applied voltage of 160 volts. Dark current values were 
higher than that stated by the manufacturer by a factor of 
mo at low applied reverse bias voltages, decreasing to a 
faerer ct 10 higher at high values. For this particular 
detector, breakdown occurs at a reverse bias voltage of 
160 volts vice the manufacturer's indicated 166 volts. 

The HP Harrison 6515A DC voltage source was then re- 
placed by a battery box consisting of a series of seven 
22.5 volt and one 1.5 volt batteries to be used with the 


avalanche detector during actual aperture averaging 
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experiments. The "Dark Current vs. Reverse Bias Voltage" 
experiment was repeated. The same equipment and oui 

as described above was used and dark current values obtained 
with the battery box were identical to those obtained with 
Exewme Harrison 6515A DC source. 

With the breakdown voltage known, the next test was to 
examine the "Signal-to-Noise Ratio vs. Reverse Bias Voltage" 
at Signal intensities similar to those previously described 
under at-sea conditions. The experimental setup was as 
Shown in Figure B3. The signal source was a .5 mW He-Ne 
(6328A) laser. The laser beam was passed through a +8 cm 
focal length lens followed by a signal chopper located just 
beyond the focal point of the lens, then proceeded a distance 
of 3.35 meters to the avalanche detector. The avalanche 
detector was mounted on a stand that permitted horizontal and 
vertical movement in order to maximize signal intensity 
appearing on the oscilloscope. The signal from the detector 
was sent to a PAR (Princeton Applied Research) Model 113 
Pre-Amp, then to a Monsanto AM-6419/USM-368 oscilloscope 
for display. External triggering was used on the oscillo- 
scope. To obtain the desired signal intensity, a Jodon 
Model PM-550 Digital Optical Power Monitor, placed at the 
location of the avalanche detector, measured signal intensity 
directly. With the uninterrupted laser beam incident on the 
Jodon detector at close range (about 10 cm), the Power 
Momitoreind cated a Maximum power output Of .39 mW. The 


diameter of the laser beam was 6 mm wide while the diameter 
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of the active detection area of the Jodon detector is 3 mm. 
Since the laser beam is theoretically (and appeared to be 
from visual inspection) Gaussian, roughly 87% of the output 
intensity was ene on the Jodon detector active area, 
covering it completely, and giving a source power output 

of 5.5X10° W/cm’. At the position of the avalanche detector, 
the intensity had decreased, due to beam expansion from the 
unfocused optics of the source output, to 1.0X107 °W/cm?. 
With the chopper added and running at normal operating 
speed, the intensity reduced further to 5.7X107‘W/cm?. 

With the diverging lens in place, to complete the setup, 
intensity decreased to 3X10° °W/em?. To reduce the intensity 
further, neutral density filters (Kodak Wrattan Filter (in 
faeeeass) N.D.4, N.D.3 and N.D.2) were used singly and in 
combination. With the N.D.4 filter inserted (reduces 
intensity by WW") between the laser and the +8 cm lens, 
Signal intensity was reduced to 3X10° W/cm“, as desired. 

When using the PAR, it was found that a saturated 
Signal occurred on the oscilloscope if PAR Gain setting 
was too high. A saturated signal was identifiable as a 
rectangular-shaped 8.75 volt signal with normally visible 
Scintillation absent. 

On 26 January 1976 an experiment comparing "S/N vs. 
Reverse Bias Voltage" at a given signal intensity for the 
avalanche detector was conducted followed by a comparison 
Oetts/N of Avalanche Detector vs. S/N of PIN-SDP Detector” 


at varying intensities. In both experiments the equipment 
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setup of Figure B3 was used. PAR settings were: L-F 

RODD OFF--100, H-F ROLL OFF--300K, Gain--Variable. The 
Harrison 6515A DC Power Supply was used to supply reverse 
bias voltage to the avalanche detector and a 6328A filter 
was used in front of the avalanche and PIN photodiode 
detectors. No lens system was used to focus the incoming 
light onto the detector active area for either detector. 
Table BI gives data collected during the "S/N vs. Reverse 
Bias Voltage" for the avalanche detector, shown graphically 
in Figure B4. In Table BI, note (1) the lower PAR Gain 


settings required to avoid saturation as V increases 


bias 
beyond 120 volts, and (2) the increase in noise voltage 


as V increases. From Figure B4 it can be seen that S/N 


bias 
increases fairly linearly from a reverse bias voltage of 
80 volts through 159 volts. Table BII is a comparison of 
the performance of the PIN photodiode detector and the 
avalanche photodiode detector. Both detectors were located 
at the same distance from the signal source and both used 
the same PAR and oscilloscope for signal output. It is 
evident that the avalanche detector, when set at its 
optimum reverse bias voltage, gives a significant improve- 
ment in performance over the PIN photodiode detector, 
especially at the weaker intensities characteristic of 
actual field experiments. 

A responsivity curve was plotted to compare experimental 
results with the manufacturer's stated performance charac- 


teristics. From Figure B5 it can be seen that experimental 
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results follow very closely the predicted values with the 
following exceptions: (1) actual Signal gains corresponding 
to manufacturer's highest predictions were not attained and 
(2) the maximum reverse bias voltage that could be applied 
was less than manufacturer indications. 

In summary, the TIXL69 avalanche detector proved to 
give a significant improvement over the PIN-5DP photo- 
diode detector for low intensity signal detection, as 
anticipated. The maximum (and optimum). reverse bias voltage 
that can be applied is 159 volts, slightly less than previous- 
ly expected. Measurements of responsivity agree with 
manufacturer's stated values but do not attain the highest 


values predicted. 
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aie bl Ca) 
Circuit Diagram: Avalanche Photodiode Detector 
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Figure Bl1(b) 
Test Setup: Dark Current Experiment 
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TABLE BI 
S/N VS. Veras 


"FOR AVALANCHE DETECTOR 


° 
26 January 1976, Source: .5mW, 6328 A He-Ne laser. 
Using chopped signal and diverging lens, as shown in 


Figure B3, signal intensity using N.D.F 4 ~ 3X107°W/cm?. 





VeTAS pon V (NOISE) V (SIGNAL) S/N 
20V 10K 55V 1.4V DAS 
50V 10K 60V TS GH 2.6 
80V 10K 65V 1.9V 2.9 
100V 10K .70V a oS 
120V 10K 1.05V 5.25V 5.0 
140V SK 1.25V 7.0V 5.6 
150V 2K 1.10V 6. SV 5.9 
155V 1K 1.00V 6. 0V 6.0 
159V 200 70V 4.35V 6.2 
160V NO SIGNAL 
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APPENDIX C 


EQUIPMENT PERFORMANCE CHARACTERISTICS 


Precision Instruments, Model No. PI-6204, four-channel 
recorder: 
Input sensitivity, .05 to 10 volts rms; 
Rueducheymenresponse, DG to 10 KHz +5 dB at 37.5 ips 
on FM Record; 


pienal to Noise Ratio, 35 dB at 37.5 ips, 300 Hz 
Oe LOUehiZ . 


Sanborn HP-3914, 14-channel Analog Tape Recorder: 
Maximum Input, 5 volts (peak to peak); 
Camee Unitey trom Oto ss0 Hz, -.25 dB at 1 KHz; 


Tape Speed 7.5 ips. 
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